AN ELEMENTARY DERIVATION OF THE MONTGOMERY 
PHASE FORMULA FOR THE EULER TOP 



JOSE NATARIO 

Abstract. We give an elementary derivation of the Montgomery phase for- 
mula for the motion of an Euler top, using only basic facts about the Euler 
equation and parallel transport on the 2-sphere (whose holonomy is seen to be 
responsible for the geometric phase). We also give an approximate geometric 
interpretation of the geometric phase for motions starting close to an unstable 
equilibrium point. 



Introduction 

The motion of the Euler top is governed by the Euler equation, whose generic or- 
bits are periodic. By applying the Stokes theorem to a suitable surface on T*S0{3), 
Montgomery |Mon91| obtained a Berry-Hannay-like formula for the angle by which 
the final position of the Euler top is rotated with respect to the initial position 
after one period (see also [MMROO, Mar92 ). A different derivation, based on the 
Poinsot description of the motion and the Gauss-Bonnet theorem, was given by 
Levi f Lev93j Fl 

The purpose of the present paper is to give a third, more elementary derivation, 
in the sense that it utilizes only basic facts about the Euler equation and parallel 
transport on the 2-sphere. The basic observation is that the motion of a fixed 
orthonormal basis as seen in the Euler top's frame can be easily understood in 
terms of the Euler flow on a sphere of fixed angular momentum norm and parallel 
transport on this sphere. 

The structure of the paper is as follows: the first section briefly reviews the theory 
of the Euler top; the main result is proved in the second section, with Montgomery's 
formula deduced as a corollary in the third section; the fourth section contains an 
approximate geometric interpretation of the geometric phase for motions starting 
close to an unstable equilibrium point. 

1. Euler top 

In this section we briefly review the theory of the Euler top, mainly to flx the 
notation. This material is standard and can be found in almost any book on 
mechanics (e.g. |Arn97l iGPSOl IMR 991 [OTi02] ) . 

An Euler top is a rigid body with a fixed point moving freely in an inertial 
frame. Its motion is described by a curve S : 5*0(3) which at each instant 

gives the orientation of the body with respect to a reference position. We have 



Partially supported by FCT (Portugal). 

^Both these derivations have been generalized to more complicated mechanical systems (see 
l AKSgSllCibOTT ). 
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S{t) = S{t)A{t), where A{t) e 5o(3) can be identified with n{t) 6 through the 
weU known hncar isomorphism defined by 

for all ^ e R'^. It is then easily shown that u:{t) = S{t)Cl{t) is the angular velocity 
of the Euler top, and hence S},{t) is the angular velocity as seen in the Euler top's 
frame. 

The (conserved) total angular momentum of the Euler top is given by 

p = s{t)m{t), 

where / : M'^ ^ is a symmetric, positive definite linear operator whose eigen- 
values /i,/2,/3 (called the principal moments of inertia) satisfy the triangular 
inequality Ii < I2 + I3 (and cyclic permutations). The 1-dimensional eigenspaces 
of / are called the principal axes of inertia. 

The vector P{t) — S~'^{t)p is the angular momentum vector as seen in the Euler 
top's frame, and is in general not constant. In fact, 

p = ^ sp + st> = ^ s{n X p) + sp ^ 0, 

and hence P{t) satisfies the so-called Euler equation 

p + r2 X p = 0. 

Using n — I^^P, this is a ODE for P{t), which in an orthonormal basis {ei, e2, es} 
of eigenvectors of / is written 

1 1 




h h 
1 1 

1 1 



P2P3 
P3P1 
P1P2 



J2 h. 

The Euler equations are integrable. An obvious first integral is ||P|| = ||p|| = p 
(which we assume to be nonzero), and a second first integral is the kinetic energy 
K=\{P,n): 

^ ^ = (P:^"'P) = (P,f^) = -{^ X P,n) = 0. 

Thus the solutions with fixed norm of the angular momentum live on the intersec- 
tions of the angular momentum sphere J|P|| = p with the ellipsoids (P, /~^P) = 2K 
for different values oi K (cf. Figure [T])0 

2. A PICTURE OF THE MOTION 

Consider a right-handed orthonormal basis {e, f , n}, fixed in the inertial frame, 
with n — p/p. Then N(^) — S~^{t)n = P{t)/p will be the unit normal vector 
to the angular momentum sphere at the point P{t), and E(t) — S'~^(t)e and 
F{t) = S~^{t) f will be (orthogonal) unit vectors tangent to the sphere at the same 
point. 



^In fact, it is easily shown that the Euler equations give the Hamiltonian flow of K on the 
sphere ||P|| = p for the symplectic structure given by — times the standard surface area element. 
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Figure 1. Integral curves of the Euler equations on the angular 
momentum sphere, assuming Ii > I2 > Is- 



We can ask ourselves what is the covariant derivative of E(i) as it moves along 
the sphere. Since e is constant, we have the Euler-like equation 

E + n X E = 0. 



As is well known |Boo03[ rdC76[ [dC93[ IO'N83| . the covariant derivative is just the 
projection E^ of E on the tangent space to the sphere. ThusH 

E^ = (E,E)E+ (E,F)F = -{ft x E, E) E - ($7 x E, F) F 

= -{fl X E, N X E> F = -(N, E X (J7 X E)) F 

= -(N, ||E|pJ7 - (E, $7) E) F = -(N, ft) F 

1 IK 
= -i(P,0)F = -— F. 

P P 

In other words, E(t) rotates with constant angular velocity —2K/p with respect 
to a parallel-transported frame on the sphere (and the same must obviously be 
true for F{t)). This gives a nice picture of the motion: the fixed orthonormal 
basis {e, f , n} is seen in the Euler top's frame as the moving orthonormal basis 
{E(i), F(t), N(t)}. If we imagine this basis placed on the point P{t) of the angular 
momentum sphere (which moves according to the Euler equation), then N(i) is the 
outward unit normal vector and E(i), F(t) rotate at a constant rate about N(t). 



^Here we use the vector identities (a, b X c) = (b, c X a) and a X (b X c) = (a, c) b — (a, b) c. 
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3. Montgomery's formula 



Except for the six equilibrium points and four heteroclinic orbits, all orbits of 
P(t) on the angular momentum sphere are periodic (cf. Figure[T]). After one period, 
the angular momentum vector as seen in the Euler top's frame has returned to the 
initial position. Therefore the Euler top has rotated by some angle A9 around p 
with respect to its initial orientation. Montgomery's formula gives this angle as 



where T is the period and A is the area of the region bounded by the orbit on the 
angular momentum sphere (with sign according to whether the orientation of the 
periodic orbit is induced by the standard orientation of this region) Q The first term, 
which is just the time integral of the component of u:(t) along p, is a dynamical 
phase. The second term, which can be thought of as a solid angle, is what is usually 
called a geometric phase. 

We can use our picture of the motion to give an elementary derivation of the 
Montgomery formula. First, it is clear that A9 exactly coincides with minus the 
angle by which E(i) has rotated with respect to its initial position after one orbit. 
Since E(i) rotates with constant angular velocity —2K/p with respect to a parallel- 
transported frame on the sphere, it will have rotated —2KT/p with respect to 
such a frame after one orbit (dynamical phase). By elementary geometry of curved 
surfaces jdC76[ [dC94|, IMor98l |Spi79| , upon returning to the initial point the parallel- 
transported frame will be rotated by an angle equal to the curvature of the sphere 
(1/p^) times the (signed) area A of the region enclosed by its path (geometric 
phase). Therefore we obtain 



which is Montgomery's formula. In particular, the geometric phase is seen as a 
consequence of parallel transport on the angular momentum sphere. 



It is possible to give an approximate geometric interpretation of the geometric 
phase in the case of an orbit which starts close to one of the unstable equilibrium 
points on the angular momentum sphere (cf. Figure[T]). For such orbits, the motion 
of the Euler top is roughly as follows: initially it is rotating with angular speed 
2K/p ~ p/h about p (that is, P(t) remains close to the equilibrium point); then 
it quickly reorients itself by rotating 180° about an axis perpendicular to p {P{t) 
rapidly moves to the other unstable equilibrium point along a heteroclinic) in 
this upside-down orientation it continues to rotate about p with the same angular 
speed 2K/p (P{t) remains close to the other equilibrium point); finally, it performs 
another reorientation by rotating 180° about a (different) axis perpendicular to p 
(P(i) returns to the first equilibrium point along another heteroclinic). In other 
words, with respect to an auxiliary frame in uniform rotation about p with angular 
speed p/h, the motion of the Euler top during one period consists basically of two 



'The two possible choices, A and —{Airp'^ — A) = A — iwp^ , lead to the same angle mod27r. 
'This is the twisting phenomenon studied in |ACC9ll ICB97) 




4. Interpretation of the geometric phase 
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180° rotations about certain axes perpendicular to p0 Since the auxiliary frame 
rotates by 2KT /p during one period, the geometric phase must result from the 
composition of the two reorientations (which in particular has to be a rotation 
about p). 

It is easily seen that the heteroclinics are half great circlesQ During the reorien- 
tations, the auxiliary frame rotates with respect to the Euler top's frame about the 
axes perpendicular to these great circles, in the direction given by the Euler flow. 
According to the geometric construction described in the appendix, the composi- 
tion of these two rotations is a rotation about the middle moment of inertia axis by 
twice the angle between the two heteroclinics. But this is exactly A/p^ , where A is 
the (signed) area enclosed by the two heteroclinics. Therefore the Euler top rotates 
by ^A/p^ about p with respect to the auxiliary frame, giving the geometric phase. 

Appendix 

Here we describe the geometric construction for the composition of rotations 
given in |Pen04j . This construction is dual to the more familiar Rodrigues-Hamilton 
construction, which uses triangles with vertices on the rotation axes |Nee97l IPar65[ 
[Sti92] . 




Figure 2. Geometric construction for the composition of rotations. 

Suppose that one wants to compose a rotation about some axis a by an angle 9 
with a rotation about some other axis b by an angle (f. Then on the unit sphere 
one draws the great circles perpendicular to a and 6, oriented according to the 
direction of the rotation. From an intersection point O of the two great circles one 
moves backwards along the first great circle by an angle |, thus reaching a point 
A, and forwards along the second great circle by an angle ^ , thus reaching a point 
B (cf. Figure [2]). Then the composition of the two rotations is a rotation about the 

^More precisely, if {g, h, n} is an orthonormal basis fixed in the auxiliary frame, so that 
g = -E- X g, then G = S'^g satisfies = - ^ ~ 0. 

''Assuming I\ > I2 > I3 they are contained in the planes P3 = i\/"r47~r/^ ^i- 
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axis perpendicular to the great circle through A and 5, from A to B, by an angle 
ip equal to twice the angle between A and B. 
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